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Abstract 

Convexity of a yield function (or phase-transformation function) and its relations to convexity of the 
corresponding yield surface (or phase-transformation surface) is essential to the invention, definition 
and comparison with experiments of new yield (or phase-transformation) criteria. This issue was 
previously addressed only under the hypothesis of smoothness of the surface, but yield surfaces with 
corners (for instance, the Hill, Tresca or Coulomb-Mohr yield criteria) are known to be of fundamental 
importance in plasticity theory. The generalization of a proposition relating convexity of the function 
and the corresponding surface to nonsmooth yield and phase-transformation surfaces is provided 
in this paper, together with the (necessary to the proof) extension of a theorem on nonsmooth 
elastic potential functions. While the former of these generalizations is crucial for yield and phase- 
transformation condition, the latter may find applications for potential energy functions describing 
phase-transforming materials, or materials with discontinuous locking in tension, or contact of a 
body with a discrete elastic/frictional support. 

Keywords: yield surfaces; phase-transforming materials; phase-transforming surfaces; nonsmooth 
of elastic energy; nonlinear elastic contact. 
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1 Introduction 



• Yield or phase-transformation functions 

Bigoni and Piccolroaz (2004) have proposed a new yield (or phase-transformation) function within the 
class of isotropic functions of the stress tensor a defined by 



in which, having defined 



Pc + C 

the meridian and deviatoric functions take the form 
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respectively, where p, q and 9 are stress invariants. 

To preserve convexity of the yield surface, the seven material parameters defining the meridian shape 
function f(p) and the deviatoric shape function g(6) are restricted to range within the following intervals 

M > 0, p c > 0, c > 0, < a < 2, m > 1, < /3 < 2, < 7 < 1. (6) 

The interest in the above yield function and in the more general class of functions {TJ lies in the 
fact that they can model the behaviour of many materials of engineering importance, such as ceramic 
(Piccolroaz et al. 2006) and metal (Bier and Hartmann, 2006; Hartmann and Bier, 2008; Heisserer et 
al. 2008) powders, metals (Hu and Wang (2005); Wierzbicki et al. 2005; Coppola and Folgarait, 2007), 
high strength alloys (for instance, Inconnel 718) (Bai and Wierzbicki, 2008), shape memory alloys (for 
instance, NiTi, NiAl, CuZnGa, or CuAINi) (Raniecki and Mroz, 2008), concrete (Babua et al. 2005), 
and geomaterials (Dal Maso et al. 2007; Descamps and Tshibangu, 2007; DorMohammadi and Khoei 
2008; Maiolino, 2005; Mortara, 2008; Sheldon et al. 2008). Moreover, eqn. (p} can be used as a general 
expression to set the condition for phase-transformations, for instance, to determine the stress threshold 
for martensitic or austenitic transformation (Raniecki and Lexcellent, 1998 and Lexcellent et al. 2002). 



• Convexity of yield or phase-transformation functions 

With reference to the class of functions (JTJ) , Bigoni and Piccolroaz (2004) have proved a general propo- 
sition providing necessary and sufficient conditions relating convexity of the yield function to convexity 
of the corresponding yield surface in the Haigh-Westergaard stress space (or principal stresses represen- 
tation) , H a crucial property in the development of new expressions for yield or phase-transformation 



1 The stress invariants p, q and 6 are defined by 



tr <r r— 1 / 3^ J 3 

P= <? = V3J2, 9 = -arccos \^— ] 

where Ji and J3 the second and third invariant of the deviatoric stress S 

J 2 = -trS 2 , J 3 = -tvS 3 , S = cr-—I, (5) 

2 3 3 w 

in which I is the identity tensor. 

2 Four years later, exactly the same proof has been independently published by Raniecki and Mroz (2008). 
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criteria. This proof is based on both the hypotheses of smoothness of the function g(9), d and validity 
of the smoothness limiting conditions g'{0) = and g'(w/3) = 0, while as noticed by Laydi and Lex- 
cellent (2009, see Appendix IA.2I for a detailed discussion), the class of functions ^ -even under the 
particularizations |(3j)- may describe deviatoric yield surfaces with corners (Fig. [TJ, in which case the 
conditions for convexity provided by Bigoni and Piccolroaz (2004) remain only necessary but not suffi- 
cient. Since the convexity proposition is fundamental in developing new yield or phase-transformation 
criteria, it has immediately attracted a strong attention (Taillard et al. 2007; Laydi and Lexcellent, 2009; 
Lavernhe-Taillard et al. 2009; Saint-Sulpice et al. 2009; Valoroso and Rosati, 2009) and may definitely 
be important in analysing yield criteria with corners. Therefore, it becomes imperative to generalize the 
convexity proposition to nonsmooth deviatoric yield surfaces, which is obtained in the present article 
(Theorem gJl Section gj. 




Fi gure 1; Examples of yield surfaces with corners obtained within the general class of yield function |£Q), namely F(cr) = 
— /c/ff( 7r /3) + l/g{&) with the Bigoni and Piccolroaz (2004) deviatoric function, eqn. (J3J2 , taking 7 = 1. (a) Deviatoric 
section: ,3 = 0,0.5,1,1.5,2. (b) Yield surface in the biaxial plane ai/ft vs. 0-2/ ft, with 0-3 = 0: /3 = 0, 0.5, 1, 1.5, 2; ft and 
f c denote tensile and compressive uniaxial yield stress, respectively. 

The generalization of the Bigoni and Piccolroaz (2004) proposition requires the generalization to 
nonregular functions of a theorem given by Hill (1968) regarding convexity of elastic strain potentials. 
In particular, Hill (1968) has shown that convexity of a smooth scalar isotropic function of a second- 
order symmetric tensor (a work-conjugate strain measure in his case) is equivalent to the convexity of the 
corresponding function of the principal values (the principal stretches in his case). The Hill's theorem is of 
fundamental importance, since in many cases [for instance for the Ogden (1982) constitutive equations for 
rubber elasticity and the so-called ' ^-deformation theory materials', Neale (1981)] constitutive equations 
of finitely-strained elastic materials are formulated with reference to the principal stretches and not with 
reference to the tensorial quantities, so that this theorem is usually reported in books (see for instance 
Ogden, 1984). Bigoni and Piccolroaz (2004) have recognized that the Hill's theorem can be useful 
also for yield functions in plasticity theory, indeed the theorem has been duplicated (with a slightly 

3 The function F{&) given by eqn. JTJ is always nonsmooth along the hydrostatic axis. However, this fact has no 
consequences on convexity, as shown in Lemma [4. II The fact that f'(p) blows up to infinity when p tends to p c and — c, 
eqn. J3ji, is the only possibility to obtain smooth closures at the hydrostatic axis. 

4 It should be noted from Fig.[T]that, although the function g(0) is smooth, the limiting conditions are g'(0) < and 
g'(ir/3) > 0, so that there are corners (yet the yield surface still results convex, see Theorem 14. l|l . 
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different proof, without mentioning Hill's theorem) in the context of elastoplasticity by Yang (1980). 
However, until now no generalization of the Hill's theorem to nonregular functions has ever been given. 
Such a generalization may be relevant for elastic strain energy functions describing phase-transformation 
materials, or for elastic potential functions describing contact with discrete elastic asperities, or materials 
with discontinuous locking in tension (Fig. [21), but it is certainly of great interest for yield functions, 
which are often nonsmooth [for instance Hill (1950), Tresca and Coulomb-Mohr] . The generalization is 
provided in Section [3] and is the basis for the subsequent generalization of the Bigoni and Piccolroaz 
(2004) proposition to yield criteria (or transforming functions) with corners (Section [4]) . 




Figure 2: Contact with a nonlinear constraint. For u < u* the behaviour is linear elastic, while at u = u* and 
u = u° stiffness becomes infinite due to contact with frictional asperities obeying a rigid-plastic Coulomb rule. Although 
the sketched model is intrinsically inelastic, a representation restricted to the loading branch could be described by a 
continuous, convex and nonsmooth potential, for which Theorem 13. f I applies. In tension, a behaviour similar to that 
sketched on the right would correspond to a discontinuous locking mechanism of a material. 



2 On smoothness of yield (or phase-transformation) functions 

The conditions for smoothness of function F{<t), eqn. |T]), can be obtained by analysing the gradient of 

d^ = -3 f{p)I+ u~ 9 s -n7 s ' (7) 

where J is the identity tensor and 

\2q' * \!3 q dcT V2Vsm30 

Note that: 

• S is discontinuous along the hydrostatic axis, but this discontinuity does not affect convexity of 
F((t), see Lemma [mi 

• S is discontinuous along the hyperplanes defined by 9 = and 9 = 7r/3. This discontinuity can 
be eliminated for functions g such that g' = at 9 = and 9 = n/3, which is the case of many yield 
functions, for instance, all yield functions described by eqn. |(3])2 when < 7 < 1. The analysis of 
the nonsmooth case, at 9 = and 9 = ir/3 is the main target of the present article and leads to 
Theorem 14.1) which is generalized into Theorem 14.21 and finally leads to Theorem 14.31 

We analyse now smoothness of the deviatoric part q/g{9) as a function of Si, S2, where Si, £2 denote 
two principal values of the deviatoric stress. Assuming that g{9) is continuous and strictly positive in 



S 2 - -q 2 I- -COS39S 
9 y 3 



(8) 
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[0,7r/3], and smooth everywhere in (0,7r/3), the gradient o{q/g{9) with respect to the variables Si,S 2 
is given by 

dq/gifi) 1 dq g'{6) 89 

1^27^7^' *=1,2, (9) 



where d 



dSi g(6)dSi g 2 (6)dS l 7 

^ = ^[^-(-DV], -1,2, (10) 

^=-^,4)™,, i = 1,2, (11) 

in which the indices are not summed and the vector m has the components: {m} = {S 2 , —Si}. 
The function H(Si,S2) is a piecewise constant function defined by 

*> " ff-'SS " " H(Sl - * )p * + * )( * + 2&)L (12) 

which takes the values 1 or —1 only, Fig. [31 We may note from Fig. [3] that the function ff(Si, S2) is 




Figure 3: Plot of the function J?, eqn. I|12|l . in the deviatoric plane. Axes a\, 0-2,0-3 are the projections of the principal 
stress axes onto the deviatoric plane. 

discontinuous along the projections of the principal stress axes on the deviatoric plane: 

axis <Ti : {Si,-Si/2,-Si/2}, if Si > then 6 = 0, if Si < then 6 = tt/3, 
axis & 2 ■ {-S2/2, S 2 , -S2/2}, if S 2 > then 6 = 0, if S 2 < then 6 = tt/3, (13) 
axis 03 : {-S 3 /2,-S 3 /2,S 3 }, if S3 > then 6 = 0, if S 3 < then 6 = tt/3. 
Accordingly, 

the function q/g{0) is smooth if and only if g{8) is smooth everywhere in (0,7r/3) and g'(0) = 
3'(tt/3) = 0, see Fig.H 



6 Note that there is a misprint in (Bigoni and Piccolroaz, 2004): their eqn. (39)i should be replaced by eqn. IllOjl . 
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3 Nonsmooth, convex and isotropic functions 

With reference to elastic potential of finite-strain constitutive equations, Hill (1968) has proven that 
convexity of a smooth scalar isotropic function of a second-order symmetric tensor (a work-conjugate 
strain measure in his case) is equivalent to the convexity of the corresponding function of the principal 
values (the principal stretches in his case). The Hill's theorem is of fundamental importance, since in 
many cases constitutive equations of finitely-strained elastic materials are formulated with reference to 
the principal stretches and not with reference to the tensorial quantities, see for instance Ogden (1984). 
Bigoni and Piccolroaz (2004) have evidenced that the Hill's theorem also applies to yield functions in 
elastoplasticity theory. 

Until now, no generalization of the Hill's theorem to nonsmooth function has ever been given. As 
mentioned in the Introduction, such a generalization may be relevant for elastic strain energy functions 
describing phase-transformation materials (although in those cases usually non-convexity is employed), 
or for potential functions describing discontinuous locking in tension, or contact with discrete elastic 
springs (as explained in Fig. [21 where contact of a rigid punch with a linear elastic set of springs having 
different heights and rigid/frictional devices is envisaged. The loading branch of this model can be 
described through a piecewise linear and convex strain energy function). 

In any case, the generalization of the Hill's theorem is certainly of great interest for yield functions, 
which are often nonsmooth, as for instance in the cases of the Hill (1950), Tresca, modified- Tresca, and 
Coulomb-Mohr yield surfaces. The generalization is provided in this Section. 

We begin with a simple lemma. 

Lemma 3.1. Let us consider a scalar isotropic function 4> of tensorial argument o~ij G Sym and the 
corresponding function <fi written with reference to the principal values Oi : 

4>( a n, °22, 0-33, o"i2, C13, 023) = </>(oi, a 2, 0-3), 
then, due to isotropy, the following equality holds 

4>(ai,a2,cr 3 ) = <f>(ai, <r 2 , 03, 0, 0, 0), (14) 
so that <p is the restriction of <fi to the subdomain of diagonal tensors. 
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Proof. The property l|14p is easily proven by the following consideration. The isotropy of 4>(<r) implies 
that the function 4>(cr) is equal to a function <f> of the invariants of cr, 

<Mctii, 022,033, 012,013,023) = <£(tr <r, tr cr 2 , tr cr 3 ) , 

and thus 

4>(o- 1 ,a 2l o- 3 ) = </>(trer,trer 2 ,tr<T 3 ) = </>(oi, o 2 , o 3 , 0, 0, 0). 

□ 

We need now to introduce the notion of subdifferential (or subgradieni) , which will be used in the 
sequel for the generalization of the Hill (1968) theorem to nonsmooth functions. 
A function <f> : U C R" — > R is convex if and only if the subdifferential 

90(X O ) = {Q e R" : <j>{X)-4>{X )>Q-{X-X ), VXeU}, (15) 

is defined and non empty at every point Xo of its domain U. 

Note that although the subgradient is a set of vectors, in the sequel we shall denote with the term 
'subgradient' both the set itself and its elements. 

The following lemma, necessary to the proof of Theorem (|3.ip , is similar to the analogous given by 
Hill (1968), but now it has been generalized and extended to nonsmooth isotropic functions. 

Lemma 3.2. Given a convex function of the principal stresses, 0(oi, 02, 03), the algebraic order of 
components of the subgradient {Qi,Q2,Qs) at (01, 02, 03) is the same as (01, 02, 03). 

Proof. From the strict convexity of <j), it follows that 

3 

J2(Q t ~Q° l )(a 1 ^^)>0, (16) 

i=l 

V(Qi,Q 2 ,Q 3 ) G 00(<7i,CT 2 ,oa) and V(Q?,Q°,Q°) e 30« a 2 °, ag). Choosing (ag,a?,^) = (01,05,03) 
and taking into account isotropy, it follows that 

(Qi-Q 2 )(<ti-o 2 )>o, (17) 

and similarly for each of the other pairs. It follows that the vector (Qi, Q2, Q3) is ordered in the same 
algebraic order as (01,02,03), a property which remains true also assuming convexity '>' instead of 
strict convexity '>'. □ 

Note also that we will make use in the proof of Theorem 13.11 of an auxiliary property of the scalar 
product between two symmetric tensors, first noticed by Hill (1968): 

"if their eigenvalues are given, but their axes are directly arbitrarily, the product attains its 
greatest value when the major and minor axes are pairwise coincident. " 

We refer to Appendix IA.ll for a detailed discussion and proof of this auxiliary property. 

Theorem 3.1. Extension of the Hill (1968) theorem to nonregular functions. Convexity of an isotropic 
(not necessarily smooth) function of a symmetric (stress) tensor cr is equivalent to convexity of the 
corresponding function of the principal (stress) values Oi (i = 1, 2, 3^. In symbols, given: 

4>(a) = 0(01,02,03), (18) 

then Vcr G Sym, 

3Q G Sym : 0(<r') - <j>(cr) > Q ■ {cr' - cr) , Vcr' g Sym, (19) 
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t 

3 

3(Qi,Q 2 ,Q 3 ) £t 3 : 0(a' 1 ,a' 2 ,a' 3 )-4>(a 1 ,a 2 ,a 3 )>Y,Qi(<r'i-<ri), V(^, 4, *£) G R 3 . (20) 

Proof. The proof that (fl9| (|20|) follows immediately from the property fl4|) . The converse (|20| =>■ 
lfl9|) is not trivial and is proven in the following. 

We denote by (01, a 2 , 03) the principal values of a given er and by (Qi,Q 2 , Q3) the subgradient of (j> 
at (cti , a 2 , (T3). We define now Q £ Sym to be 

Q = QiQi ®q 1 + Q 2 q 2 ® <72 + <23<? 3 ® <?3: 

where q 2 j 93} is an orthonormal basis of R 3 . Then, assuming that (a[, a' 2 , cr' 3 ) are numbered in the 
same algebraic order as (cti, 02, 03), and since from the Lemma [3721 we know that the algebraic order 
of {Qi,Q 2 ,Q3) is also the same as (cti, 02, 0-3), the auxiliary property of the scalar product (proven in 
Appendix I A. ip . implies that 

3 

^2 Qi - Vi)>Q- (<?' - o-), Vcr' G Sym. (21) 

i=l 

Since, by hypothesis, the following equation holds true 

3 

<t>{tr') - <j>(<r) = a' 2 , a' 3 ) - fa, a 2 , cr 3 ) > ^ - o*) , VK,^,^) G R 3 , (22) 

i=i 

eqn. f2Tj) guarantees that Q G d4>{cr), so that 0(c) results to be convex. □ 

4 Convexity of yield functions with corners 

We begin with proving that the discontinuity of the yield function gradient along the hydrostatic axis 
[see eqn. JH])] is inconsequential on convexity. 

Lemma 4.1. The convexity of the function q/g(0) is unaffected by the fact that S and S defined in 
eqn. (E) are discontinuous along the hydrostatic axis, where Si = S 2 = 0. 

Proof. From the definition of convexity, it follows that (van Tiel, 1984) 

q/ g{9) is convex at (S±, S 2 ) for every line t through (Si, S 2 ), the restriction to t 

of q/g{0) is convex 

Let us consider all deviatoric lines through the point {Si = 0,S 2 = 0}, these can be represented 
(using parameter e and slope k) as {e,ke}. The restriction of q/g(9) to these lines is a function h(e), 
whose derivative with respect to e is 

h'(e) = l^--q^e\ -{l,k}, (23) 

where V denotes the gradient taken with respect to the variables {Si, S 2 }, so that, since V0 • {1, k} = 
and 

q=\e\V3y/l + k + k*, Vg = ^ S ; lgn ' ={2 + fc,l + 2fc}, (24) 
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we obtain 

\/2> sign e 



h'(e) = V y Vl + fc + fc 2 . (25) 
At a singular point, convexity requires that 

lim h'(e) < lim h'(e), (26) 

which is always satisfied, so that the discontinuities in S and S along the hydrostatic axis are incon- 
sequential on convexity. □ 

With reference to the deviatoric part q/g{0) of the yield function ((J), we give now necessary and 
sufficient conditions (Theorem 14. 1|) for equivalence between convexity of yield functions and convexity 
of yield surface. To this purpose, we first need the following lemma. 

Lemma 4.2. Given a generic isotropic function (j> of the stress that can be expressed as 

$(a 1 ,a 2 ,a 3 )=$(S l ,S 2 ), (27) 
where Si and S 2 are two of the principal components of the deviatoric stress, i.e. 

Si = - (2d! - a 2 - a 3 ) , S 2 = ^ {-<t x + 2a 2 - 0-3) , (28) 
convexity of </>(o"i, 02, 03) is equivalent to convexity of 4>{S\,S 2 ). 

Proof. This proposition follows immediately from the fact that the relation ([28)1 between {Si,S 2 } and 

{ci, (T2, 03} is linear. □ 

The following theorem is the generalization of Lemma 3 by Bigoni and Piccolroaz (2004) to the case 
of nonsmooth deviatoric sections of the yield surface. Note that the difference between the two versions 
of the theorem lies on the two conditions g'(0) < and (/(7r/3) > 0. A consequence of the following 
theorem is that the convexity conditions by Laydi and Lexcellent (2009) are only sufficient (but not 
necessary) for convexity of smooth functions (see Appendix IA.2|) . 

Theorem 4.1. Convexity of nonsmooth deviatoric representation q/g(9) vs. convexity of the deviatoric 
section of the yield surface. 

Assuming that g{9) is continuous and strictly positive in [0,7r/3] and twice- differ entiable everywhere 
in (0,7r/3), convexity of 

W) (29) 

as a function of S\,S 2 is equivalent to the convexity of the deviatoric section in the Haigh-Westergaard 
space: 

g 2 + 2g' 2 - gg" > 0, V6»e(0,7r/3) and g'(0) < 0, g'(7r/3) > 0. (30) 

Proof. Let us define ft as the set of all points {Si,.^} not on the axes cti, <t 2 ,(T3, see eqns. (fl3|) . The 
theorem is proven first (point 1 below) by showing local convexity at all points of Q (regular points) 
and, second (point 2 below), considering the points of dft, i.e. the axes &\ , a 2 , 0-3 , where the function 
has corners (singular points). 

1) Local convexity in f2 (for which < 9 < ir/3). 

The function q/g{9) is C 2 (0), so that we can apply the convexity criterion based on the Hessian. 
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The Hessian of the function (l29l) is 



d 2 g/g(9) _ 1 

dSidSj <r 



d 2 q fn , 2 „. 96» 961 , / a 9 89 dq 39 d 2 9 



q(2g' -gg") aa ao - gg' flc ac + ac ao +g 



^ v a 'dSidSj aa \dStdSj dSjdSi ^dSidSj 

' (31) 

where i and j range between 1 and 2 and all functions q and are to be understood as functions 
of Si and S2 only. The Hessian of q may be easily calculated to be 



d 2 q 27 



dS.dSj 



where indices are not summed and {m} = {S2, — Si}. The Hessian of 9 becomes 
d 2 9 -1 / cos 36* d cos 36* d cos 39 d 2 cos 39 s 



dSidSj 3 sin 36» V sin 2 36* dS t dSj dSidS 3 

so that 

dq 39 dq 39 d 2 9 -1 



dSidSj dSjdSi H dS t dS J 3 sin 36* 
wher^l 



d 2 qcos39 nn d 2 q cos 39 d cos 39 d cos 39 
cos 30 



dSidSj dSidSj %in 2 30 dSi 8S 3 

(32) 



<9cos36> 9-\/3sin30- <9 2 gcos36> 2 J 3 

H{Si,S 2 )mi, = -27 —m.mj. (33) 



dS^ 2q 2 v " z/ " SSidSj 

A substitution of ((33)) into ((32]) yields 

<9g 89 dq 89 8 2 9 
dSldS] + dS~8S~ + q dS l dS J 

so that we may conclude that the Hessian l(3T]l can be written as 



= 0, (34) 



d 2 g/g(9) 27 (g 2 + 2g' 2 -gg") 

^ds~ = - 7? mimj ' (35) 

from which condition g 2 + 2g' 2 — gg" > is immediately obtained. 

2) Local convexity on <9£1 (for which 9 = or 9 = tt/3). 

We consider in the following only the axis <ti , since the proof remains strictly similar for the other 
axes. 

2.1) Case = 0. 

A line t through (1, —1/2) has the parametric representation {(1 + e, — l/2 + fce)|e G R}, where 
e is the parameter and k is the slope of the line. Using this representation, the restriction to 
t of q/ g{9) is a function h(e), whose derivative is given by eqn. (|23| . From the limits 

1 1 ok 

q^3/2, Vq^ {-3/2,0} and Vfl -> ± {1,2} as e -> ± , (36) 

|1 + 2k\ v 3 

we derive 

Ai( „,^ T | |1+2t| v|. (37) 

so that the convexity condition for h(e) at e = 0, namely, h'_(0) < h' + (0) is equivalent to 
g'(0) < 0. 



Note that there is a misprint in (Bigoni and Piccolroaz, 2004): their eqns. (43)i should be replaced by eqn. (1331) 1 . 
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2.2) Case = tt/3. 

A line t through (—1, 1/2) has the representation {(—1 + e, 1/2 + ke)\e e R}, where k is the 
slope of the line. Using this representation, the restriction to t of q/g(9) is a function h(e), 
whose derivative is given by eqn. (j23"l) . From the limits 



1 + 2k 

q^3/2, Vq^ {-3/2,0} and V6 -> T — — -— {1,2} as e^O*, (38) 
we derive 

i4(<»~£±^|i + u|$j|, (39) 

so that the convexity condition for /i(e) at e = 0, h'_(0) < h' + (0) is equivalent to g'{i:/3) > 0. 

Since q/g{9) is locally convex in the sets fl and 912, the proof is concluded by noting that ft U dQ 
represents the whole deviatoric plane, so that q/g{9) is globally convex. 

□ 

There are yield surfaces, for instance that proposed by Hill (1950), see Fig. 0(a), presenting corners 
for values of internal to the interval (0, 7r/3). In particular, assuming a piecewise smooth function 
g{9), the Hill (1950) criterion can be formulated within the general class of yield functions {TJ, namelj0, 
introducing the yield stress in uniaxial compression — f c , by F{cr) = —f c + q/g(9), where 



cos 



cos 



icos 1 (cos 30) 



7T 1 

cos- 1 (cos30) 

3 3 



< 9 < tt/6, 

(40) 

tt/6 < 9 < 



Another example of a deviatoric section with corners in 9 = 0, 9 = tt/3, and 9 = 9\ = 77r/30 is given 
by F{&) = -fc/ 9 (n/3) + q/g(0), with 



cos [tt/12 - l/3cos _1 (cos30 1 )] 



< 9 < 



cos [tt/12 - 1/3 cos- 1 (cos 30)1 

9(fi) = { r , , (41) 

cos [tt/4 — 1/3 cos 1 (cos30i)J 



cos [tt/4- 1/3 cos- 1 (cos 30)] ' 



h < 9 < n/3, 



which is plotted in Fig. [5] (b). 

It is clear from the above examples that employing the function g{9) defined by eqn. with 
different values of parameter f3 on a finite number of subintervals of G [0, tt/3] , it is possible to represent 
all possible nonsmooth deviatoric sections of a yield surface. This statement justifies the interest in the 
following theorem, covering the situations in which the yield surface presents corners for values of 
internal to (0,7r/3). 

Theorem 4.2. Convexity of piecewise- smooth deviatoric representation g(9) vs. convexity of the deviatoric 
section of the yield surface. 

Assuming that g{9) is continuous and strictly positive in [0,7r/3] and twice- differ entiahle almost ev- 
erywhere in (0,7r/3), and denoting by 9i € (0,7r/3) the singular points of g(9), convexity of 

(42) 



7 Bigoni and Piccolroaz (2004) have noted that the Hill criterion cannot be expressed by the function g(S), eqn. (J3J2 , 
defined on the whole interval 9 £ [0, 7r/3] through a unique value of parameter (i. 
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Figure 5: Yield surface deviatoric sections presenting corners for values of 9 internal to the interval (0, 7r/3): (a) Yield 
criterion proposed by Hill (1950), described by eqn. ||40|I . (b) Yield criterion described by eqn. JHJ. 



to the convexity of the deviatoric section in the Haigh-Westergaard 

(43) 

(44) 

Proof. Conditions lf43j) have been already proven in Theorem 14. II and do not need further explanation. 
We therefore restrict our attention to the singular points 9i, to derive condition (|44|) . 

We consider a generic point in the first 7r/3-sector of the deviatoric plane (taken clockwise from axis 
<7i; the proof can be easily extended to the other sectors), 2q/3{cos9, — cos(7r/3 — #)}, and the parametric 
representation (with the parameter e) of all lines of slope k through this point (Fig. [6]) 



as a function of Si, S2 is 
space: 

g 2 + 2g' 2 - gg" > 0, V# G (0, tt/3) - {<?.;} and g'(0) < 0, g'{ir/3) > 

and 



-q cos ( 



e,--gcos 



(45) 



The derivative of the restriction h of q/g{&) to this line is again given by eqn. I|23p . with all the functions 
calculated at points (|45|l . Taking the limit values at e = 0, 



J3 „ V3 . 
Vg = < — cos a — sir 



, -V3sin6»|, V6» = —h\cos(^ - 6>) , cos 6>| 



and noting that H = — 1 in the 7r/3-sector under consideration (see Fig. [3]), we obtain 
V3 



lim h'(e) - 



V3 cos 9- (l + 2fc) sin 9 



where 



± for —00 < k < 



2 g*(6) 

cos (7r/3 — 9) 



(l + 2fc)cos6» + V3sin6l 



(*) = 



=F for 



cos (7r/3 — 9) 



cos( 



cob 9 



< k < +00. 



(46) 



(47) 



(48) 
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Using eqn. (|47|l into condition h'_(0) < h' + (0) yields in both cases inequality lj44"|). 



□ 




Fi gure 6: Bundle of lines {cos 6 + e, — cos(7r/3 — 9) + ke\ in the deviatoric plane. Note that k represents the slope of the 
lines in the nonorthogonal reference system 6"i,5"2. 



We are now in a position to state the generalization of the Proposition 1 given by Bigoni and 
Piccolroaz (2004) to yield surfaces with corners. 

Theorem 4.3. Convexity of piecewise-smooth yield function vs. convexity of the yield surface. 

Convexity of the yield function (lj) is equivalent to convexity of the meridian and deviatoric sections 
of the corresponding yield surface in the Haigh- Westergaard representation. In symbols: 

r r>o, 

g 2 + 2g' 2 - gg" > 0, V# € (0, tt/3) - {9 t } 
<?'(0)<0, g'(ir/3)>0, 



convexity of F(tr) = f(p) 



9(0) 



(49) 



where g(9) is a continuous and strictly positive function in [0,7r/3], twice- differentiate almost everywhere 
in (0, 7r/3), and Qi denotes the singular points of g(6). 



Proof. The proof follows directly from Lemma 14.21 and Theorems 13.11 and 14.2 



□ 



5 Conclusions 

Yield surfaces used in elastoplasticity theory often have corners. For these nonsmooth functions, we have 
given in this paper a general theorem providing necessary and sufficient conditions for the equivalence 
between the convexity of the deviatoric yield function and its representation as a surface in the Haigh- 
Westergaard stress space. This theorem is useful for the definition of new yield function or transformation 
function for phase-transforming materials. We have also provided a generalization to nonsmoothness of 
a theorem relating convexity of a scalar isotropic function of tensorial variable to the convexity of the 
corresponding functions of the tensor principal values. This can find applications in the formulation of 
nonsmooth-convex elastic potential energy functions. 
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A APPENDIX 



A.l Proof of the auxiliary property of the scalar product of two symmetric 
tensors 

We provide the proof of the auxiliary property of the scalar product of two symmetric tensors, which is 
often used (among others, by Ogden, 1984). The property has been noticed by Hill (1968), who did not 
provide a complete proof (which is only sketched in a footnote), perhaps because of a lack of space. We 
were not able to find a proof of the property anywhere. 

Theorem A.l. Let A, B be two symmetric tensors. Then, denoting by cti,a 2 ,a 3 and Pi, {3%, fa the 
eigenvalues of A and B, respectively. 

A ■ B < aipx + a 2 (3 2 + a 3 p 3 , (A.l) 

given that the eigenvalues of the two tensors are numbered in the same algebraic order. 

Proof. Given the eigenvalues of the two tensors, we keep the eigenvectors a±, a 2 , 03 of A fixed and seek 
for the maximum of A- B as the eigenvectors bi 7 b 2 ,b 3 of B rotate with respect to ai, a 2 , a 3 . Therefore, 
the problem can be formulated in terms of the following optimization problem 

max A ■ B, (A. 2) 

b 1 ,b 2 ,b 3 

with the constraint that (61, 62, 63) be an orthonormal basis, 

bhi ■ b N = S MN , (A. 3) 

where Smn is the Kronecker symbol. 

This optimization problem can be solved using Lagrangean multipliers, so that we maximize the 
function 

A B = ai/3i(ai ■ bi) 2 + ai[3 2 {a x • b 2 ) 2 + aif3 3 {a x ■ b 3 ) 2 
+a 2 /?i(a 2 • bi) 2 + a 2 (3 2 (a 2 • b 2 ) 2 + a 2 /3 3 (a 2 • b 3 ) 2 

+a 3( 3i(a 3 • bi) 2 + a 3 /3 2 (a 3 • b 2 ) 2 + a 3 /3 3 (a 3 ■ b 3 ) 2 (A.4) 
+Ai(6i ■ 61 - 1) + A 2 (b 2 • 6 2 - 1) + A 3 (6 3 • 63 - 1) 
+A 4 (6i • b 2 ) + A 5 (b 2 ■ 63) + A 6 (6 3 • bi), 
as a function of bi,b 2 , b 3 and the Lagrangean multipliers Aj (i = 1, 6), thus obtaining 

dA ■ B 



dbx 
dA-B 

db 2 
dA ■ B 

db 3 



2/3iA6i + 2Ai6i + A 4 6 2 + A 6 6 3 = 0, 
2f3 2 Ab 2 + 2A 2 b 2 + A 4 bi + A 5 6 3 = 0, (A.5) 
2f3 3 Ab 3 + 2A363 + A 5 6 2 + A 6 b! = 0, 



together with the constraints (|A.3|1 . 

In the case of distinct eigenvalues 0i , (3 2 , (3 3 , the system (| A.5[) is satisified if and only if 

6 M • Ab N = 0, for M ^ N, (A.6) 



1C 



and thus if and only if 61, 62, &3 are eigenvectors of A. This proves that the extreme values of A ■ B are 
attained when the two tensors are coaxial. The maximum is then selected from six possibilities. 

In the case Pi = fh ^ 03, the same line of thought used above allows us to conclude that the extreme 
values of A ■ B are attained when 63 is an eigenvector of A, in which case the two tensors A and B are 
coaxial and, choosing 63 = a 3 , A B = (ai + 02)^1 + a 3 (3 3 . 

The case /3i = {3% = (3$ is trivial. The two tensors A and B are coaxial and the scalar product is 
A ■ B = (qi + a 2 + a 3 )Pi. □ 

A. 2 The convexity condition given by Laydi and Lexcellent (2009) 

Laydi and Lexcellent (2009) have shown, with the example reported in Fig. [TJ that the convexity con- 
ditions given by Bigoni and Piccolroaz (2004) does not cover yield surfaces with corners. In fact, by 
selecting within the class (U) the following function 

F(a) = -f t + q/g(9), g~\e) = 2 - cos 2 6, (A.7) 

the Bigoni and Piccolroaz (2004) conditions are satisfied, but, although g'(0) = 0, the yield surface has 
concave corners at 6 = 7r/3, g'(ir/3) < [instead of g'(ir/3) > 0, corresponding to convex corners], see 
Fig. [3 




Fi gure 7; The example by Laydi and Lexcellent (2009) showing the possibility of describing corners within the class of 
functions {TJ with the choice l|A.7ll . Since g' (71-/3) < 0, hypotheses of Theorem 14.11 are violated and indeed the deviatoric 
section of the yield surface has reentrant corners at 9 = tt/3. (a) Deviatoric section, (b) Yield surface in the biaxial plane 
cri /ft vs. 1T2/ ft, with 0-3 = 0. 

Laydi and Lexcellent (2009) incorrectly argued that the Bigoni and Piccolroaz (2004) proposition 
on convexity had flaws, while the problem lies only in the fact that the deviatoric section of the yield 
surface described by eqns. (|A.7|I has corners, a case which is not covered by the Bigoni and Piccolroaz 
(2004) proposition and has been addressed in the present paper. 

Laydi and Lexcellent (2009) also provided sufficient conditions for convexity of the deviatoric section 
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of a smooth yield surface. These conditions in our notation read 

- cos 9(gg') + sm9g 2 > 0, 

cos9(gg') + sm0(2g' 2 - gg") > 0, 



(A. 



for all 6 G [0,tt/3]. 

However, these conditions are neither necessary, nor sufficient for deviatoric sections with corners, 
while the correct, necessary and sufficient conditions are those specified by Theorem 14. 11 To fully justify 
this statement, we provide the two counter-examples below. 

Counter-example 1: Conditions l|A.8j) are not necessary for convexity of yield functions, even with 
smooth deviatoric section. 

This is made clear by the following counter-example (taken from eqn. ^2 with (3 = 0.5 and 7 = 0.99): 



F(a) 



fc 



g(w/3) g{0y g{6) 



0.5 



7T cos" 1 (0.99 cos 36») 
6 3 



(A.9) 



The deviatoric shape function l|A.9j) 9 corresponds to a smooth and convex deviatoric section, see Fig. 
[8j but it is easy to show that it does not satisfy the condition (|A.8Jl o. 



g(K/3) = 
(smooth) 



(a) 



g'(O) = 
(smooth) 




3/2 o 3 // c 



(b) 





2 /f t 1-5 








up. 








/ 0.5 






-1.5 f 


-1.0 -0.5 


0.5 


y.O 1.5 




-0.5 








-1.0 








-1.5 







Figure 8: An example of convex and smooth yield surface, corresponding to eqns. |)A.9|| . not satisfying the Laydi and 
Lexcellent (2009) conditions l)A.8ll . (a) Deviatoric section, (b) Yield surface in the biaxial plane <ri/ft vs. 02/ ft, with 
oz = 0. 



Counter-example 2: Conditions (|A.8(1 are not sufficient for convexity of deviatoric sections with 
corners. 

This is made clear by the following counter-example: 



F(a) 



fc 



9(0) 



9(0) = e 2 -o. 



- V sin ( 



(A.10) 
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g(n/3)<0 
(non-convex comer) 



rn 



(a) 



g'(0) = 
(smooth) 




3/2 o 3 // B 




Figure 9: An example of non-convex yield surface, corresponding to eqns. llA.lOjl . satisfying the Laydi and Lexcellent 
(2009) conditions llA.811 . (a) Deviatoric section, (b) Yield surface in the biaxial plane a\j jt vs. o"2//t, with 0-3 = 0. 



The deviatoric shape function l|A.10|) 9 corresponds to a non-convex deviatoric section, see Fig. 02 but 
it is easy to show that it does satisfy both conditions l|A.8|) . 
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